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High-energy behavior of amplitudes in a gauge theory can be reformulated in terms of the evolution 
of Wilson- line operators. In the leading logarithmic approximation it is given by the conformally 
invariant BK equation for the evolution of color dipoles. In QCD, the next-to-leading order BK 
equation has both conformal and non-conformal parts, the latter providing the running of the 
, coupling constant. To separate the conformally invariant effects from the running-coupling effects, 

f****) ■ we calculate the NLO evolution of the color dipoles in the conformal A/"=4 SYM theory. We define the 

"composite dipole operators" with the rapidity cutoff preserving conformal invariance. The resulting 
£SJ , Mobius invariant kernel for these operators agrees with the forward NLO BFKL calculation of Ref. 
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The high-energy scattering in a gauge theory can be described in terms of Wilson lines - infinite gauge factors 
Qh| ordered along the straight lines (see e.g. the review Q). Indeed, the fast particle moves along its straight-line classical 
trajectory and the only quantum effect is the eikonal phase factor acquired along this propagation path. In QCD, 
for fast quark or gluon scattering off some target, this eikonal phase factor is a Wilson line - an infinite gauge link 
ordered along the straight line collinear to particle's velocity n M : 
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I. INTRODUCTION 



U v (x±) = PexpUg [ du A^{un + x±)\, (1) 

J — oo 



Here is the gluon field of the target, x± is the transverse position of the particle which remains unchanged 
throughout the collision, and the index r\ labels the rapidity of the particle. 

The high-energy behavior of QCD amplitudes can be studied in the framework of the evolution of color dipoles. 
Let us consider the small-x behavior of structure functions of deep inelastic scattering (DIS). At high energies the 
virtual photon decomposes into quark and antiquark which propagate along the straight lines separated by transverse 
distance and form a color dipole - two- Wilson-line operator. 



U"(x ± ,y ± ) = 1 - ±i X {fr>(x ± )fr'>{v ± )} (2) 

The energy dependence of the structure function is then translated into the dependence of the color dipole on the 
rapidity 77. There are two ways to restrict the rapidity of Wilson lines: one can consider Wilson lines with the 
support line collinear to the velocity of the fast-moving particle or one can take the light-like Wilson line and cut the 
rapidity integrals "by hand" . While the former method appears to be more natural, it is technically simpler to get 
the conformal results with the latter method of "rigid cutoff" in the longitudinal direction. 

Thus, the small-x behavior of the structure functions is governed by the rapidity evolution of color dipoles 0,0]. 
At relatively high energies and for sufficiently small dipoles we can use the leading logarithmic approximation (LLA) 
where a s « 1, a s \tixb ~ 1 and get the non- linear BK evolution equation for the color dipoles 



^-W( Zl ,z 2 ) = jd 2 z 3 ^^[W{z 1 ,z 3 )+U\z 3 ,z 2 ))-U\z ll z 3 )-W(z 1 ,z 3 )W{z 3 ,z 2 )] (3) 

Z 13 Z 23 
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where r) = ln^p- and z i2 = zi — zi etc. (As usual, we denote operators by "hat"). The first three terms correspond 
to the linear BFKL evolution Q and describe the parton emission while the last term is responsible for the parton 
annihilation. For sufficiently low xb the parton emission balances the parton annihilation so the partons reach the 
state of saturation [8] with the characteristic transverse momentum Q s growing with energy f /xb (for & review, see 

It is easy to see that the BK equation ([3]) is conformally invariant in the two-dimensional space. This follows from 
the conformal invariance of the light-like Wilson lines. Indeed, the Wilson line 

U(xx) = Pexp Ug [ dx + A + {x + ,x ± )\ (4) 

J — oo 

is invariant under the inversion x^ — > x^ /x 2 (with respect to the point with zero (-) component). Indeed, (x + , x±) 2 — 
—x\ so after the inversion x± — > x\_jx\ and x + — > x + jx\ and therefore 

r f°° x + x + -\ 

U(x±) -> Pexp Ug d— A+( — ,x±) \ = U{x ± /x 2 ± ) (5) 

It is easy to check that the Wilson line operators lie in the standard representation of the conformal Mobius group 
SL(2,C) with conformal spin (see appendix A). 

The NLO evolution of color dipolc in QCD [Kj] is not expected to be Mobius invariant due to the conformal anomaly 
leading to dimensional transmutation and running coupling constant. However, the NLO BK equation in QCD [101 ] 
has an additional term violating Mobius invariance and not related to the conformal anomaly. To understand the 
relation between the high-energy behavior of amplitudes and Mobius invariance of Wilson lines, it is instructive to 
consider the conformally invariant J\f = 4 super Yang-Mils theory. This theory was intensively studied in recent 
years due to the fact that at large coupling constants it is dual to the IIB string theory in the AdS$ background. In 
the light-cone limit, the contribution of scalar operators to Maldacena- Wilson line [ll[ vanishes so one has the usual 
Wilson line constructed from gauge fields and therefore the LLA evolution equation for color dipoles in the Af = 4 
SYM has the same form as ©. At the NLO level, the contributions from gluino and scalar loops enter the picture. 

As we mentioned above, formally the light-like Wilson lines are Mobius invariant. Unfortunately, the light-like 
Wilson lines are divergent in the longitudinal direction and moreover, it is exactly the evolution equation with respect 
to this longitudinal cutoff which governs the high-energy behavior of amplitudes. At present, it is not known how to 
find the conformally invariant cutoff in the longitudinal direction. When we use the non-invariant cutoff we expect, 
as usual, the invariance to hold in the leading order but to be violated in higher orders in perturbation theory. In our 
calculation we restrict the longitudinal momentum of the gluons composing Wilson lines, and with this non-invariant 
cutoff the NLO evolution equation in QCD has extra non-conformal parts not related to the running of coupling 
constant. Similarly, there will be non-conformal parts coming from the longitudinal cutoff of Wilson lines in the 
Af = 4 SYM equation. We will demonstrate below that it is possible to construct the "composite conformal dipole 
operator" (order by order in perturbation theory) which mimics the conformal cutoff in the longitudinal direction so 
the corresponding evolution equation has no extra non-conformal parts. This is similar to the construction of the 
composite renormalized local operator in the case when the UV cutoff does not respect the symmetries of the bare 
operator - in this case the symmetry of the UV-regularized operator is preserved order by order in perturbation theory 
by subtraction of the symmetry- restoring counterterms. 

Let us present our result for the NLO evolution of the color dipolc in the adjoint representation (hereafter we use 
notations Zij = z% — Zj and {T a )b c = —if abc ) 

^0imr nl (6) 
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where 
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[Tr{l>^}] conf = Tr{UlU%} + / d 2 z 3 -|^[Tr{r"^^T»^ 3 ^ } _ JV c Tl{C%t)£}] In (7) 



2tt 2 
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is the "composite dipole" with the conformal longitudinal cutoff in the next-to-leading order and a is an arbitrary 
dimensional constant. (Similar expression for the conformal two-dipole operator in the r.h.s. of this equation is 
presented below, see Eq. ([48])). In fact, a(rj) = ae v plays the same role for the rapidity evolution as for the usual 
DGLAP evolution: the derivative -4- gives the evolution equation ©. The kernel in the r.h.s. of Eq. © is obviously 

2 2 2 2 

Mobius invariant since it depends on two four-point conformal ratios z %> z % 4, and z \ 2Z l i . We will also demonstrate that 

r—, Z 14 Z 2? 1 Z 13 Z 24 

Eq. ([6|) agrees with forward NLO BFKL calculation of Ref. |l[ . 

The paper is organized as follows. In Sect. [IT] we remind the derivation of the BK equation in the leading order 
in a s . In Sect. IIII1 which is central to the paper, we calculate the scalar and gluino contributions to the small- x 
evolution of color dipoles and assemble the NLO BK kernel in TV = 4 SYM. We rewrite the NLO BK kernel in the 
conformal form © in Sect. IIVI and compare our results with the NLO BFKL calculations in Af = 4 SYM in Sect. 
fVl . In Sect I VII we derive the NLO BK equation for the composite dipoles in the fundamental representation, both 
in JV = 4 SYM and in QCD. The Mobius group for the light-like Wilson lines is presented in Appendix A while in 
Appendix B we find the explicit form ([7]) of the composite conformal dipole operator by calculating the appropriate 
impact factor. In Appendix C, we find the leading-order evolution for the four- Wilson-line operator and Appendices 
D and E contain some technical calculations which may distract readers from main discussion. 



II. DERIVATION OF THE BK EQUATION 



Before discussing the small-x evolution of color dipole in the next-to-leading approximation it is instructive to recall 
the derivation of the leading-order (BK) evolution equation. 

For the NLO calculation we use the lightcone gauge p^A^ — 0. In addition we find it convenient to use the "rigid 
cutoff" prescription in the longitudinal direction. We consider the light-like dipoles (in the p\ direction) and impose 
the cutoff on the maximal a emitted by any gluon from the Wilson lines so 

/>oo 

U% = Pexp ig dup± A^upi + x±) 

K( x ) = I -\a k \)e- lk - x A,{k) (8) 



(hereafter we use the ^.-inspired notation ct n p = m^jz for brevity) . Note that the cutoff © respects the unitarity of 
Wilson lines {U^U^ = 1). 

2 

The momenta p\ and p2 are the light-like vectors such that q — p\ — xbPi and p — p2 + — pi where p is the 
momentum of the target and m is the mass. Throughout the paper, we use the Sudakov variables p — api + (3p2 +p± 
and the notations x m = x^p^ and = x^p 1 ^ related to the light-cone coordinates: x* = x + x, = x~ yj s/2. 

Our metric is (1,-1,-1,-1) so p 2 — a(3s — p\ and x 2 — -x.x* — x\. 

To find the evolution of the color dipole ([2]) with respect to the slope of the Wilson lines in the leading log 
approximation we consider the matrix element of the color dipole between (arbitrary) target states and integrate 
over the gluons with rapidities r\\ > rj > 772 = 771 — A77 leaving the gluons with 77 < 772 as a background field (to be 
integrated over later). In the frame of gluons with r\ *~ 771 the fields with 77 < 772 shrink to a pancake and we obtain 
the four diagrams shown in Fig. [1] Technically, to find the kernel in the leading-ordrer approximation we write down 
the general form of the operator equation for the evolution of the color dipole 

±Tr{UlU%} = K hQ Tr{UlU%} + ... (9) 

(where dots stand for the higher orders of the expansion) and calculate the l.h.s. of Eq. (J9j> in the shock-wave 
background 

shockwavc — shockwavc 

(10) 

In what follows we replace (...)shockwavc by (...) for brevity. 

With future NLO computation in view, we will perform the leading-order calculation in the lightcone gauge p^A^ = 
0. The gluon propagator in a shock-wave external field has the formpjl, [IH 
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FIG. 1: Leading-order diagrams for the small- a; evolution of color dipole. Gauge links are denoted by dotted lines. 
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where 



dp„(k) = fid, - —(k^p 2 u + k^p 2fJl ) - —p2fj.P2u 



(12) 



Hereafter we use Schwinger's notations (x±_\F(jp±)\y±) = Jdp e l ( p,x ~ v > J -F(px) (the scalar product of the four- 
dimensional vectors in our notations is x • y = -(x*y, + x*y,) — (x,y)±). Note that the interaction with the shock 
wave does not change the a-component of the gluon momentum. 
We obtain 



g 2 du dv{At^{up 1 +x 1 )A b ^{vpx+y A _)) Fi& m = -Aa s / — ( x± \^L^U ab ^^\y ± ) (13) 



J -co 



Formally, the integral over a diverges at the lower limit, but since we integrate over the rapidities rj > 772 we get 
(in the LLA) 



g 2 I du I dv (Ar n (up 1 +x ± )A b . m (vPi+y±)) Fi M = -^sA V (x ± \^U ab ^\y ± ) 
' > P±_ P± 



and therefore 



(ui^u^M = ~^v(t a u Zl ®t b ui)[d 2 Z3 { -?^u« b 



(14) 



(15) 



13^23 



The contribution of the diagram in Fig. [lb is obtained from Eq. £5) by the replacement t a U Zl ®t b U\ 2 -> U Zl t b (g)UU, a , 
z 2 Z\ and the two remaining diagrams are obtained from Eq. [14] by taking z 2 = %i (Fig. \Tfc) and z\ — z 2 (Fig. [lji). 
Finally, one obtains 



« I — I 7T ./ 2.1 o 2no 



13^23 



-13 



,2 ^23 

z 23 



(16) 



so 



{^{uim)l\ g m =^p-[ d 2 z 3 ^mT«u zl ulxu z M 2 } - ^{u zi ut 2 }} (17) 

61_1 7T J z 13 z 23 iV c 

There are also contributions coming from the diagrams shown in Fig. [2] (plus graphs obtained by reflection with 
respect to the shock wave) . 
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FIG. 2: Leading-order diagrams proportional to the original dipole. 

These diagrams are proportional to the original dipole Tr{U Zl Ul 2 } and therefore the corresponding term can be 
derived from the contribution of Fig. 1 graphs using the requirement that the r.h.s. of the evolution equation should 
vanish in the absence of the shock wave (when U = 1). It is easy to see that this requirement leads to 

{Tr{U£Ufr}) = ^ fd 2 z 3 ^[TT{T a U Zl UlT a U Z3 Ut 2 }-N c TT{U Zl Ul 2 }} 

~ J z 13 z 23 

which gives the BK equation for the evolution of the color dipole in the adjoint representation:. 

±Tr0lU^} = J / d 2 z 3 -^[Tr{T a U Zl UlT a U Z3 U} 2 } ~ N c Tr{U Zl Ul 2 }] (18) 
arj it j z is z 23 

III. CALCULATION OF THE NLO BK KERNEL IN Af = 4 SYM 

In the next-to-leading order the contributions to the kernel come from the one-loop diagrams for the color dipole 
in the shock- wave background. We will take the results for the gluon part of the NLO BK kernel from Ref. [10] and 
calculate the contribution of scalar and gluino loops. We use the Af = 4 Lagrangian in the form (see e.g. Ref. [14]]): 

C = - -F^Fpu + -(.D"*?) - -g 2 f abc / ,,BC *f*^$^*™ 

+ KA°fV^f - * K A ^A B A s b \^ k f abc + i K A X sAB Af\% c r bc (19) 

Here $j are scalars, X" A gluinos and EJ 7 = (t]ab>^Vab)j ^ij = (Vab>~ ^Vab) where rj^ B are standard 't Hooft 
symbols. The bare propagators are 

/p—ip-(x — y) p j^y) qV 1 

2M . , (>tf(x) \V(y)) = / <r^HP-(x-») J*£a (20) 
p z + ie H J p z + ie 

and the vertex of gluon emission in the momentum space is proportional to (k\ — k2)^T a 8u for the scalars and <r M T a 
for gluinos. (We do not need Yukawa or four-scalar vertices at this level). The diagrams in the shock- wave background 
are calculated similarly to the tree diagrams discussed in the previous Section. 

A. Gluon contribution to NLO BK 

Let us start with the gluon contribution to the NLO evolution kernel. There is no difference between the gluon 
part of the kernel in QCD and in Af = 4 SYM so we will just copy it from Ref. [Hi) replacing tx{t a U Zl t b U\ } in the 
fundamental representation by Tr{T a U Zl T b W Z2 } (th roughout the paper we denote traces in the fundamental and the 
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adjoint representations of color group by tr{...} and Tr{...}, respectively) 
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where [i is the normalization point in the MS scheme. Our normalization fro Gell-Mann matrices is ti{t a t b } = \8 ab ■ 
Note that the last term in r.h.s. is Mobius invariant. The coefficient stands in front of the non-conformal 
terms coming from the running of the coupling constant and as discussed in the Introduction, there is an additional 

non-conformal term ~ In t 2 - In -j 2 - coming from the non-invariance of the longitudinal cutoff (TSl) . 

z 13 Z 23 

It should be noted also that there is one small difference between QCD and N = 4 calculations of the gluon loop 
due to the fact that in supersymmetic theories it is more natural to use the dimensional reduction scheme [15j | instead 
of dimensional regularization. In dimensional reduction scheme the factor g^g^" = d± coming from the product of 
three-gluon vertices should be replaced by 2: g^ v g^^ v — > 2. Making proper replacement in formulas in Sect. IV of 
Ref. [l(| one gets the factor ^ in the r.h.s. of the above equation, in place of ^ in Eq. (5) in Ref. [l(|). 



B. Contribution of scalar particles 

1. Diagrams with two scalar- Shockwave intersections 
First, we calculate the diagram with two scalar-shockwave intersections shown in Fig. [3J 




FIG. 3: Diagrams with the scalar loop bisected by the shock wave. 
The scalar propagator in the shock- wave background has the form Q : 



($?(*)$$(!/)) = e(x*y*)6 al >- ctact[3(x ± \ 



i(a/3s — pj_ + ie) 



\y±) 



da 
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^)^(-y*)(xx|e- i ^ x «^ ab e i ^^|y J .)+^(- a; *)&(^)(a ;± |e^^c7t a6 e-^^|yx) 



We start with the calculation of the Fig. [3^ diagram. Multiplying two propagators (I22p . two scalar-gluon vertices 
and two bare gluon propagators we get 
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(USU^m = g 2 / d« / dw(i?(«pi + zi)iS(«pi + «a)) (23) 
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where u = a\/a and it = 1 — it. The contribution of the diagram in Fig. is obtained by replacing e^ 9l+92,Zl ^ by 
_ e *(9i+<?2,z2K anc j the two remaining diagrams in Fig. [3); and Fig. [3ji are obtained by change x <-> y. We get 
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Following the method suggested in Refs. [Hi Oil E3] we separate the UV-divergent part by adding and subtracting 
Zi -> z 3 contribution: U™' Uf' = (U™'U™' - z 4 -> z 3 ) + U™'U%f. We get 
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(30) 



The second (UV-divergent) part should be calculated at d_L ^ 2. As in the case of gluon loop, the Fourier transform 
(|2T|) at d± ^ 2 is complicated so it is convenient to return to Eq. (|26[) in the momentum representation. After 



(31) 



replacing {/™ m {J™ by [/™ m £/™ 3 ™ , integrating over u and changing variables to &2 =92 = fc' 5 P = Qi + <Z2, ' = Qi — ki 
(so that qi = p — fc', ki = p — I — k' and ki + = p — I) the Eq. (f26|) turns into 

(Tr{^f7i 2 "}) Fig [3] Z4 _3 
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(p-k'- 1) 2 - fc' 2 



fc' 2 

(p,p-0 /2 , Z 2 8 



r 2 (i-|)r(§) 



2np 2 (p-lf 



4vr(3 - e)r(2 - e) p 2 (p-l) 2 
\np 2 /l 2 \n(p-l) 2 /l 2 



(P 



2-e 



\P~1\ 



2-£ ;2-e 



/ 2 " £ ) 



3(p-0 2 



3p 2 



0(e) 



Subtracting the pole in e corresponding to counterterm (see the discussion in Refs [ToL fl^|) we obtain 
a±(Tr{UZU%}) Figm ^ Z3 = 3^N C J d 2 z Tr{T a U Zl U} 3 T b U Z3 Ui} J <t*-*p i 2 ^l 



(32) 



(33) 



2(p,p-0 



-3j> 2 (p-/) 2 

Using the Fouriers integral from Appendix A to fiol ] we get 



, I 2 8 
In- +3 



lnp 2 /Z 2 \n(p-l) 2 /f 



d 2 pd 2 l e ^(p.*is)-<(p-J,a 8 ) 
21n 



2(p,p-0 / , / 2 8 



-In- 



3p 2 (p-0 2 V p 2 



3(P~0 2 
lnp 2 /7 2 ln(p-Z) 2 // 2 



3p 2 



(213,^23) foi„ Z 13 Z 23P 2 , 16 



'13*23 



127T 2 Z 2 . 



13*23 



In 



*12 

2 2 2 
Z 13 Z 23^ 



-12 



12tt 2 ; 



23 



Z 12 



3( P -l) 2 
1 



3p 2 



127T 2 2 2 3 



In z 2 3 fi 2 + 



z 2 

— \n^f- 

z 13 z 12 
1 



12tt 2 ; 



ln z 2 3l u 2 



23 



and therefore 



iqs. 

-(Tt{UIU^}) sc ^. s 



2tt 3 



d z 3 



'13*23 



In 



2 2 2 
*13*23M 



'12 



T±{T*U z JJlT b U«pt J 



(34) 



(35) 



Combining Eqs. (|30|) and ([35]) we obtain the full contribution of diagrams in Fig. [3] to the NLO kernel in the form 

d_ 

drj ' 



3a 2 f d 2 z 3 d 2 Z4 L 
4^4 



-31 



'13*24 



Zi aZ- 



14*23 



2 ~2 



13*24 



7 2 2 
^1 4*' 



In 



14*23 



z 13 z 24 
2 ,2 



Z,aZ 



14*23 



(EC' C>i 6 ' - *4 -> z 3 )Tr{[T a 1 T b ]Ul [T a , T° ]U Z2 } 



2tt 3 



13*23 



In 



222 

Z 13 Z 23A* 



-12 



Tt{T a t%&%T a uM) 



(36) 



where we have promoted the shock-wave Wilson lines to operators. 
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2. Scalar loop 



Besides diagrams with the scalar loop bisected by the shock wave calculated above, there are diagrams with the 
ordinary scalar loop shown in Fig. 2J 





FIG. 4: Diagrams with bare scalar loop. 



The integral for the scalar loop has the form 



cf 4 ~ £ fc' 



(2k' -k)^2k> ~k) v 
ik' 2 + ie)[(k - k') 2 + it 



(4n) 2 ~iT(4-e)(-k 2 )i [ " " 9 ^> 



48tt 2 



2,18 



(k^k u -5^, v ) (37) 



and therefore the contribution of the diagram shown in Fig. |4] takes the form 



±(Tr{UlU^}) Figm = _ d 2 z 3 [Tr{T^U zl UtT a U Z3 Ut 2 } ~ ^-Tr{f/ Zl C/tj 

2 



x fct 2 M 2 k'd 2 q [e^ z ^ - e i(9.*«i)][ e -i(Mi3) _ e -i(MM)] %-§(ln ^ + ~i 

/ ' K 2 (7 2 I K 2 3 J 



(38) 



As usual we should add diagrams obtained by the reflection of diagrams shown in Fig. 3] with respect to the shock- 
wave line. Their contribution is obtained from Eq. (|38p by the replacement q «-> k in the logarithm so the final result 
for the sum of all diagrams of Fig. 2] type has the form 



(39) 



x [ct 2 kd 2 k'd 2 q [e^-) _ e i(«,*»)ir e -i(Mi3) _ e *(*«)]^(]n-^ + -) 
7 1 Arg 2 L <7 2 k 2 3 J 



16 



r z 2 z 2 
In 13 f - 


16" 


+ 


■ i 


1 " 






~ z 13 


Z 23 J 



z 23 J 
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The total contribution of scalar particles to the NLO kernel from the diagrams in Fig. [3] and [4] is a sum of Eqs. (|39|) 
and l[5g|) 



3V 
4tt 4 



d 2 z 3 d 2 



z 13 z 24 



'14^23 



hi 



'13*24 



"14^23 



: '13^'24 

2 2 
Z 14 Z 23 



(40) 



47T 3 



d 2 z 3 Tr{r a (7 2l {/t 3 r a (7,3[/t 2 }{ 



z 13 z 23 



In z 12 ^' 



z 4 - ^{[^T^UJJT^T 6 ^} 



In 



°23 ■ 



Finally one needs to add the contribution of diagrams without scalar-shockwave intersection shown in Fig|5l They 
are proportional to the "parent dipole" Tr{U zl Ul 2 }, and their contribution can be found from Eq. (|4"0")) using the 
requirement that the r.h.s. of the evolution equation must vanish as z\ — > zi (since lim Zl ^ Z2 U z 1 1 U r] \ Z2 = 1, see the 
discussion below Eq. (J8j) ) . It is easy to see that the following formula for the total contribution of scalar particle 
fulfils this requirement: 




FIG. 5: Diagrams without scalar-shockwave intersection. 



(Tr{^^ })scala 



3a 2 f d 2 z 3 d 2 Zi 
I^ 1 



°:S4 



2 _|_ Z 13 Z 24 + z 14 z 23 j n Z 13 Z 2A 



2 ^2 ^y2 ~ 2 ^2 -^2 

z 13 z 24 z 14 z 23 z 14 z 23 



(41) 



x {U? 3 aa ' Ul'f -z 4 ^ z 3 )Tr{[T a 7 T b }Ul[T a ' 7 T b '}U^} 



eg.,, 



z 3 [Tv{T a UlU^T a U Z3 Utn ~ N c Tr{UlU%}]{^ 



13^23 





8" 




- 1 


1 " 


In z 2 2 p, 2 - 


3- 


+ 




~ T 2 " 






- z 13 


Z 23 J 



In 



6 13 

z 23 



Note that we have written this equation in the operator form by promoting the shock-wave fields in the r.h.s. of Eq. 
(|40f to operators. 



C. Gluino contribution 

The diagrams for the gluino contribution to the NLO kernel are shown in Fig. [5] The gluino propagator in the 
shock-wave background has the form 



(T\% (x) \¥ (y)) = 6(x*y*) / d 4 k " a & 5 ab S u 

k + le 



-ia(x—y). 



+ 5v9(x*)6(-y*) da , 
Jo 2 " s 

roc e ia(x-y). 

+ 6ij6(-x*)9{y*) / da 
Jo 



(a:±|(opi + p ± )e- tI ^ x -'p 2 U ab e l! ^ v *{ap 1 +p±)\y±) a c 



2a 2 s 



(x x \(api +p±)e i! ^ x 'p2U^ ab e i! ^ y '(api + p±)\y±) a a 
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FIG. 6: Diagrams with the gluino loop. 



where p adl = p^ a6l and p aa = p^cr^ a . 

This propagator has the same form as the quark propagator in the shock-wave background in QCD so one can use 
the result for the quark part of the NLO BK kernel in QCD calculated in Refs. Eg. We replace tr{t a U Z3 t b UlJ 
for quarks in the fundamental representation by 2Tr{T a C/ 23 T f> [/j 4 } for gluinos in the adjoint representation (and 
tr{t a U Zl t b Ul 2 } by Tr{T a U Zl T b Ut 2 } for Wilson loop as usual) and obtain 



drj 



luino 



2ajN c 
3tt 3 



d 2 z 3 [T±{T a U^U%T a UlU%} - N c Tr{UlU%}] 



13*23 



■ ( In z\ 2 n 2 



-23 



'13*23 



Z 23 



2a 2 f d 2 z 3 d 2 Z4 



°34 



^14 2 '23 + Z ?4 Z 1 



'24*13 



2 2 2 
z 12 z 34 j n ^13*2 1 



(42) 



2(^13^24 



Z li Z 2i) 



'14*23 



}(Ul) aa '(Ul - Ul) bb 'TT{[T\T b ]Ul[T a \T b ']U^} 



D. The N=4 kernel 



Now we are in a position to assemble the NLO BK kernel in TV = 4 SYM. Adding the gluon contribution (f2Tj) . the 
scalar part (|4"Tj) . and the gluino term (|4"2")) we obtain 



dr\ 



d z 3 



'13*23 



l-^rS+21n#ln4 



4?r 



-13 



z 



23 



8tt 4 



-34 



12*34 



'13*24 



'23*14 



hi 



7 2 2 
*13*24 



'14*23 



[Tr{T a C/^C/^ r «^ 3 ^ } _ N c Tr{UlU%}) 

(43) 



z 3 <-> z 4 



x TT{[T\T b }UlT a 'T b 'u^ +T b T a Ul[T b \T a )(j^}(&l - U^) bb ' 



z 12 z 34 , 



'12*34 



Zi oZn 



111 



'13*24 



'14*23 



z 3 



2 4 



Tr{ [T a , T b ] ti\ T a 'T b ' U% + T b T a [T b ' . T a '] U% } (0£ ) bb ' } 
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Using Eq. (|114|) from the Appendix E one can rewrite this equation as follows: 



J-Tr^C/t,} 



n- 

4tt 4 



d 2 z 3 



2 r 

§r{i 



a s JV c 



2 

Z 13 Z 23 



d 2 Z3<i 2 Z4 Z 2 2 Z 2 4 r 



4tt 
1 + 



-34 



7T 

T 



7 2 2 

21n^ln^ 



2 



12 



-12 



} [Tr{T"UlU%T a UlU%} - ^Trj^^}] 



7 2 2 
z 12 z 34 



, z 13 z 24 
2 2 
Z 14 Z 23 



x Tr{[T a , T b ]UlT a 'T b 'u^ + T b T a Ul [T b \ T a ']U%}{Ul) aa ' 01 - Ul) bb ' (44) 



All terms in the r.h.s. of this equation are Mobius invariant except the double-log term proportional to In -t 2 - In t 2 -. 

z 12 z 12 

As we discussed in the Introduction, the reason for this non-invariance is the cutoff in the longitudinal direction which 
violates the formal invariance of the non-cut Wilson lines. 

It is worth noting that conformal and non-conformal terms come from graphs with different topology: the conformal 
terms come from 1— >3 dipoles diagrams (see Fig. 6 in Ref. [Iol |) which describe the dipole creation while the non- 
conformal double-log term comes from the 1— >2 dipole transitions (Fig. 9 in Ref. [lQ()) which can be regarded as a 
combination of dipole creation and dipole recombination. 



IV. CONFORMAL DIPOLE AND CONFORMAL NLO KERNEL 



A we discussed in the Introduction, it is possible to define the composite conformal dipole operator order by order 
in perturbation theory in such a way that the evolution equation for this operator would be Mobius invariant. The 
form of this operator can be guessed from the expression (|9"9"|) (see the discussion in the Appendix B) 

[Tr{l>^}] C ° nf = Tr^t^^/^^ (45) 

Z ~ J Z Yi Z Ti Z 13 Z 23 

Let us find the NLO evolution kernel for this operator and demonstrate that it is conformal. 
For the evolution of the composite operator (|4"5"j) we get 

= ^{IT^I + ll^ (46) 
Writing down the evolution of the four- Wilson-line operator (|105| calculated in Appendix C we obtain 



d 



.,2 



-12 



TT 2 J Z3 Z\ Z Z\ 3 r ' 47T 
.,2 



7T 2 Z 2 Z 2 

ZL_21n^ln^ 

3 z 12 Z 12 



} [Tr{T a UlU%T a UlU%} - N c Tr{UlU%}} 



^ / d 2 z 3 d 2 z 4 {Ul) aa ' \^l 2 T Ti{T a T b UlT a 'T b 'u^ + T b T a U , J i T b ' T a ' U^}(2Ul - IJ\ - U\) bb ' 



,2 
'13 



Z i aZ 



■Ti{T a T b Ul[T al ',T b ']U% + [T b ,T a ]UlT b 'T a 'u^}(2Ul ~ - U^) bb ' 



14^34 

-^h r TT{[T a 1 T b }U r J i T a 'T b 'u^ + T b T a Ul[T b ' ,T a ']U^}(2Ul - - U!> 3 ) bb ' 

'-2 1 '3 1 

a 2 s N c 



12 In 12 

2 2 2 2 

2 13 Z 23 z 13 2: 23 



, ^^(^{r 8 ^^^^}-^^^}) In + ... (47) 

J z 13 z 23 z 14 z 24 z 13 z 23 



where the dots stand for the last (conformal) term in the r.h.s. of Eq. (|44|) . 

Next we need the "counterterms" converting the four- Wilson- line operator Th:{T a U% Ul v T a U% UJ V } into the confor- 
mal operator. In principle, this should be done similarly to obtaining the "conformal dipole" ([7| in Appendix C: one 
should expand the T-product of conformal operators in the next (a 2 ) order in perturbation theory and rearrange the 



13 



6- Wilson-line operators in such a way that the coefficient in front of the combination [Tr{T a U^Ul^T a U^ 3 U^}] coni is 
conformal. Since it means the calculation of the NNLO impact factor which is a formidable task, we will use another 
method to get the four- Wilson-line conformal operator. We will make a guess 

{Tr{T a UlU^T a U r J 3 U^} - N c Tr{UlU^}] cont 
= Tr{T a UlU^T a UlU^} - N c Tv{U^U^ 2 } 

r 2 2 

+ FT / d 2 z i (UD aa '\TT{T a T b UlT a 'T b 'u^ + T b T a U^T b ' T a ' 'U%}(2U% - U\ - \J\) bb ' -^r In (-^r) 

8lT J I ~ ' z 14 z 24 z 14 z 24 

2 2 

_ T T { T a T b Ul[T a ' ,T b ']U% + [T b ,T a ]Ur>T b 'T a 'u%}(2Ul - (j^ - tj^) bb ' l n (-f^f-) 

z 14 z 34 z 14 z 34 
2 2 

- Tr{[T a ,T b ]UlT a 'T b 'u^ +T b T a Ul[T b \T a )(j^}(2Ul - - U^f 3 ' ln (-^§-) } 

z 24 z 34 z 24 z 34 ^ 

_ ^/^^(It{r»CT, 1 C^2«C^^J-JV c M^ 1 Ujj)ln^ (48) 

Z7T J z 14 z 24 Z 14 Z 24 

and check that it leads to the conformal evolution equation ©. 

Rewriting the evolution equation (|47[) in terms of conformal operators (|45[) and (|48p we obtain 



rfr/ 



[Tr{^ lL >I 2 "}] 
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= —n / a Z 3 



6 12 



-13^-23 
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47T 
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_ ZL _ 21n^ln^ 
z m 2 m 2 

° z 13 z 23 



+ 
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8tt J z 13 z 23 ^ 



66' ^12 



^2 »2 



+ 4(J7£) ao '(J7£ - U!! 3 ) bb ' Tr{ [T a , T b ] C/^ T a ' T b ' U}^ + T b T a [T b ' ,T a '] } 



-"23 l / ^14^23 A 
r 2 _2 V r 2 r 2 / 
z 24 z 34 z 12 z 34 
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z 24 z 34 



12*34 
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"13 



'14*34 
v2 



ln( 



-"13*23 ' 
^2 -,2 ~2 / 
z 12 z 14 z 34 



(#? a )°°'(#? a + UD bb 'Tr{[T a , T b ]UlT a 'T b 'u^ + T b T a U^ [T b ' ,T a ']U^}^_ ln ( 



'24*34 



*23*13 ' 
r 2 r 2 r 2 ; 
z 12 z 24 z 34 



J z 13 z 23 z 14 z 24 z 13 z 23 



(49) 



Note that with our accuracy we do not need to specify the form of "counterterms" for the conformal composite 
operators in the a 2 term in the r.h.s. of this equation. 

Using Eq. (|107[) and the integral (|1 16[) from Appendix we get 



di] 



conf 
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"s / j2 
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"13*23 
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13 



-23 



} [Tr{T a UlU%T a UlU%} N c Tr{UlU%}Y 



^ / d 2 z 3 d 2 z 4 2 Z 2 2 2 0l) aa '(Ul - Ul) hb 'Ti{[T a , T b ]UlT a ' T b '0% + T b T a iJ\ [T b ' ,T a ']U^} ln -^§4 

7 z 13 z 24 z 34 " z 12 z 34 ' 

a s N, 



2tt 3 



d 2 z 3 In ln - ^2C(3)Tr{^ i C/t 2 >7 } 
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-13 -°23 
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Z 13 Z 24 Z 34 



{Ul) aa '{Ul - UD bb 'Tr{[T a ,T b }UlT a ' T b ' + T b T a U^[T b ' ,T a ']U^}\n- 



14^23 
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We see that the non-conformal term double-log term ~ In In -j 2 is canceled with the correction coming from the 

2 13 Z 23 

subtitution of the dipole by the composite operator ([7|). This confirms our expression (|7|) for the conformal dipole 
and justifies our guess for the conformal composite 4- Wilson-line operator (|48p. 

Substituting the dots in the r.h.s. of this equation for the last (conformal) term in Eq. ((44)) we get the final 
evolution equation for the conformal composite operator cited in the Introduction ([5]): 

±[Tr{UlU%}] c ° nf (51) 



a. / ,, z 
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2 
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a,N r it 2 
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2 / o 2 2 

7 Z 13 Z 23 
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a s J J2„ J2„ Z 12 J 2 | n Z 12 Z 34 



[Tr{T a UlU%T a U Z3 U%} - Ar c Tr{^^ } ] 



conf 



72 72 Z 12 / 

a z^d z^— — ^ — 21 
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7 2 2 
z 12 z 34 



In 



i 13 /i 24 



"14^-23 ' 



4tt 4 

x Tr{ [T a , T b ] T a ' T b ' + T b T a [T b ' , T a '] } [[U% 3 ) aa ' (U% A ) bb ' - (z 4 -» z 3 )] 



13^24^34 
l ' T b ' JJ^ 

At this point we would like to discuss the gauge invariance of our evolution equation (|5"Tj) . The Wilson- line operator 
is gauge invariant up to gauge rotations at ±00. As it was discussed in Refs. 0,0], the evolution equation should be 
reformulated in terms of gauge-invariant Wilson loops. In particular, Tr{C/J U^} in the l.h.s. of this equation should 
be promoted to 

Tr{U^U^} = Tr{[oop 1 + z 1 ,-oop 1 + z 1 }[z 1 ,z 2 ]-oo[-oop 1 + z 2 ,cop 1 + z 2 ][z 2 ,z 1 } 00 } (52) 

where we use the notation [z%, z 2 ]± QO = [z\ ± oopi, z 2 ± oop\] and the precise form of contours connecting these points 
does not matter since the fields at infinity are pure gauges. We do not have a simple way to introduce these gauge 
links at infinity to the r.h.s. of Eq. (|51[) in the adjoint representation, but it can be easily done if one rewrites 
the adjoint traces in terms of traces in the fundamental representation using Eq. (|106p from the Appendix D. For 
example, tr{U z JJl 1 }tv{U z . 2 Ul 3 U Zl Ul 2 U Z3 Ul i \ should be replaced by 

Tr{[oopi + z 4 , -oopi + z 4 ][z 4 , zi]_ 0O [-c»pi + Zi,oopi + Zi][zi, z^} 

x Tr{[oopi + z 2 , -oopi + z 2 ][z2,z 3 ]-oo}[-oopi + z 3 ,oopi + z 3 ][z 3 , z^oolzi, z 2 ]-oo[-oopi + z 2 ,cop 1 + zi\ 

X [z 2l Z 3 ] oo [OOPI + Z 3 ,~OOpi +Z 3 ][z 3 ,Z 4 ]_ 00 [-OOpi +Z 4 ,OOpi + Z 4 ][z 4 , Z 2 ]oo} 

and similarly for other traces in the r.h.s. of Eq. (|107[) . With this replacement, the evolution equation (f5T|) is gauge 
invariant. 

V. COMPARISON TO NLO BFKL 

In this section we compare our kernel with the forward NLO BFKL results for AA = 4 SYM [l[. To compare to 
the BFKL amplitude of gluon-gluon scattering at high energies we need to expand our Wilson lines up to two-gluon 
accuracy. We define the analog of Eq. (1) in the adjoint approximation: 

1 ... N 2 - * N 2 

V&V) = 1 ~Jp—[ Tr { U SU i y r '}= j^^[U(x,y)+U(y,x)-U(x,y)U(x,y)} ~ J ^—[U( x , y )+U(y,x)} (53) 

The corresponding conformal dipole operator in the BFKL approximation has the form 

*0 = V"(zuz 2 ) + ^ A^-Ar ln ^r\V\z u z 3 ) + ^ (z 2 ,z 3 )-V"(z 1; z 2 )\ (54) 

47r J z 13 z 23 z 13 z 23 

Using color traces (|110j) . (|111[) . and (| 1 09[) from Appendix D it is possible to demonstrate that the conformal 4- Wilson- 
line operator (|48f reduces to the sum of three conformal dipoles: 

[MT^tft^^ = -^(^-lJ^^^aJ+^^a.^-Vlrf^i,^)] (55) 



and therefore the evolution equation (|51|) turns into 

jj^Lf ( Zl , z 2 ) =^J d 2 z 3 [l - K' onf (z 1; z 3 ) + V^a, z 3 ) ~ Vl lf (z 1; z 2 )] (56) 
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It is convenient to change Z4 <-+ z 3 in the second term in square brackets and to perform the integral over z& for the 
second and the third terms. One obtains: 
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(57) 



as it is easily seen from Eqs. (|114[) and Q125P from Appendix E. 

Using Eq. (|57[) one can rewrite the linearized equation (|56[) in the following way: 
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(58) 



In this form the evolution equation (158")) coincides with the conformal NLO BFKL equation restored in Ref. [HI from 
the eigenvalues calculated in Ref. |l|. 

VI. EVOLUTION EQUATION IN THE FUNDAMENTAL REPRESENTATION 

A. In AT = 4 SYM 

For comparison with QCD let us calculate the evolution equation for color dipoles in the fundamental representation. 
The gluon part can be taken from Ref. [n| 
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The scalar part can be obtained from Eq. (|4Tj) by the replacement T a — > t a and Tr 
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Similarly, for gluino contribution one gets from Eq. ((42 
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Adding together Eq. (|59[). (f60|) and (161|) we obtain the evolution equation for the color dipole with longitudinqal 
cutoff (|HJ) in the fundamental representation: 
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due to Eq. (|114p from Appendix D. 

The composite conformal dipole in the fundamental representation can be obtained from Eq. ([7])) by the usual 
substitution T a -> t a and Tr -> tr: 
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Similarly, the conformal 4- Wilson-line operator 
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Repeating the steps which lead us to the Eq. (|5Tj) in Sect. IIVI we obtain the conformal evolution equation in the 
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fundamental representation 
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Note that it can be obtained from the equation in the adjoint representation ([51]) by same replacement T a 
-> tr. 
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B. In QCD 



It is instructive also to present the evolution equation for composite operator in QCD. The resulting equation 
will not be Mobius invariant because of the running coupling constant so composite operators (|64[) and (|65p are not 
strictly speaking conformal. We will, however, keep the notation [...] conf as a reminder that these operators were 
conformal in TV = 4 SYM. 

To get the evolution equation for "conformal" composite operators (|64| in QCD we subtract the scalar l[60|) and 
gluino (|6"Tj) contributions from Eq. (j6"B"|) and add the quark contribution calculated in Refs. [HI, [l7|. We obtain 
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where b = ^-N c — \nf and we have ^ instead of ^ because in QCD we use dimensional regularization rather than 
dimensional reduction scheme. Following the analysis of Ref. [lOj we will outline how the above kernel reproduces 
the NLO BFKL eigenvalues 

In the two-gluon approximation the conformal dipoles (|64[) reduces to 
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and therefore the evolution equation (|67[) turns into 
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where we used formula 
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following from integrals (|114|l and (|125|l from Appendix E. 

For the case of forward scattering (li(x,y)) = U{x — y) and the linearized equation ([70| can be reduced to an 
integral equation with respect to one variable z = z\i- Using integrals (104)-(106) from Ref. [l0( and the integral 
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The function — fa~rp In 2 + -PX?> + $(<Z, <?') enters the NLO BFKL equation in the momentum space [l9[ and 
since the eigenfunctions of the forward BFKL equation are powers both in the coordinate and momentum space, it is 
clear that the corresponding eigenvalues coincide. As to the first term in the r.h.s. of Eq. I|72[). one can demonstrate 
using the analysis carried out in Ref. [10( that this term also agrees with the eigenvalues of Ref. |l9j. However, this 
analysis would lead us away from the main topic of this paper so we defer it until our next publicaton. It should 
be also mentioned that the statement in our previous paper [l(| that our equation disagrees with NLO BFKL 
was due to erroneous calculation of the integral (|1 14|) which was assumed to be zero. After taking into account the 
(5-function contributions in the r.h.s. of eq. (jl 14[) the disagreement disappears. 



VII. CONCLUSION 



The amplitudes in N = 4 SYM are conformally invariant and therefore the Regge limit (fSTj) of these amplitudes 
must be invariant with respect to Mobius transformations of the transverse plane. If we want to use the operator 
expansion to find this amplitude, it is better to expand in operators which are Mobius invariant. As we demonstrate 
in Appendix A, the light-like Wilson lines are formally invariant. However, they are divergent in the longitudinal 
direction, and at present the regularization of this rapidity divergence which respects the conformal invariance is 
not known. We manage to circumvent this problem by using the non-invariant "rigid cutoff" ([8} and restoring the 
conformal invariance order by order in perturbation theory by subtracting the proper counterterms (made again of 
Wilson lines). The resulting NLO evolution equation for "composite conformal dipoles" is Mobius invariant and agrees 
with forward NLO BFKL calculations in Af = 4 SYM. 

Let us comment on the non-conformal result of the calculation of NLO BFKL kernel in Af = 4 SYM carried 
out in Ref. [2(|. We think that the difference between our kernel and that of Ref. [2(| is due to different cutoffs 
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for longitudinal integrations. The authors of Ref. [2CJ] propose that the transformation of their kernel of the type 
-?^nlo ~^ A'nlo — &s[Klo, O] with some suitable operator O will restore conformal invariance. This is exactly what 
happens in our case of the kernel (|44|) with the "rigid cutoff" ([8|) of the rapidity divergence. Let us discuss the 

transformation proposed in Ref. [201 ] in our language. If we define V{z) as 

V n (z) = V v (z) + a s fd 2 z'f(z,z')V^z') (73) 



we see that the (linear) evolution equations for the operators V v and V v are related by 

-^nlo = -Knlo - a s [K w ,f] (74) 
Indeed, differentiating the operator (|73[) with respect to r\ we obtain 

^(,) = Av,(, ) + as / d v/(z ) ,o|vV) 

r dz'K LO (z, z')V^(z') + a s fdz'K NLO (z, z')V^(z') + a s J dz'dz" f(z, z')K h0 (z', z")V\z') 
r dz'K LO {z,z')V r '(z) 

+ a s I dz'K NJj0 (z, z')V V) + a s [dz'dz"[f(z, z')K h0 (z', z") - K LO (z, z')F{z> ', z")]V^(z) + 0(a 3 s ) (75) 



which corresponds to Eq. (|71)l . Our transition between V v and V! 7 f is of the type (|73[) so it is not surprising that 
the kernels (|44|) and ([6]) are different. We think that one can recover the conformal kernel ([6|) from the kernel of Ref. 
pol ] as long as one finds the appropriate /(z, z'). It should be also mentioned that the transformation (|73|) with both 
V 1 ' and f(z, z 1 ) conformally invariant does not change -Knlo as can be easily seen from Eq. ([75]) . Thus, the form of 
the conformal composite dipole is not unique (our Eq. ([7]) is one particular choice) but the conformal kernel -/^nlo is 
unique. 

In conclusion let us discuss possible generalizations of our method. The operator expansion of the type ()82jl is 
relevant for processes like deep inelastic scattering where the strong gluon fields come from the nucleon (or nucleus) 
target and the spectator (virtual photon) is a weak source of the gluon field. For the processes like heavy-ion collisions, 
the projectile-target symmetric language of 2+1 - dimensional effective action seems more adequate (here 2 is the 
number of transverse dimensions and 1 stands for rapidity). There are many attempts in the literature to find 
comprehensive effective action [2l| but the answer for the ultimate high-energy effective action eludes us so far. It 
is possible that considering this problem in TV = 4 SYM where we have the additional requirement of conformal 
(Mobius) symmetry to restrict the effective action will help us to find the correct effective action at high energies. 
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VIII. APPENDIX 
A. Conformal properties of the light-like Wilson lines 

In this Section we demonstrate that the light-like Wilson lines are invariant under the conformal (Mobius) group 
SL(2,C). It is easy to demonstrate that the set of transverse-space operators 

'S-^'-i^-tK 2 ), ~S^\{D-iM 12 ), § + ^ l -(P 1 + l P 2 ) (76) 
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form an SL(2,C) algebra: 

[o0)&t] = ±S±, [S+, S-] = 2Sq, [Sq,S±]=±S±, [S+, S-] = 2S (77) 

Here we use standard textbook definitions of the momentum operator P, angular mpmentum operator M, dilatation 
operator D, and special conformal generator K . Using the conventional commutators of these generators with the 
gluon field 

i[P", A a ] = d^A a , i[D, A a ] = {x^ + l)i Q , 
i[M" v , A a ] = (x^d u - x v d»)A a - {g ua A^ - g^ a A v ) 

i[Ki*,A a ] = (2x fi x u d" - x 2 d» + 2x»)A a -2x v {g va A> i - g» a A v ) (78) 
one can easily obtain the action of these operators on the light-like Wilson lines. In complex notations 
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: : : ■■ = ■■■-'■■'-■ ~ -(— '—,)■ ol 2{JU ' (79) 



1 , „.J2 - 1 „-„2 

dz 2' 



these commutators take the form 

[S-, U{z, z)] = z 2 d z U(z, z), [So, U(z, z)] = zd z U(z, z), [§+, U(z, z)] = -d z U(z, z) 
[L, U(z, z)] = z 2 d z ll{z, z), [So, U(z, z)] = zd- z U(z, z), [§+, U(z, z)\ = -d z U(z, z) 

These equations mean that the operators U (z, z) lie in the standard representation of conformal group SL(2, C) with 
the conformal spin 0. 

B. NLO impact factor 

As we demonstrated in the Appendix A the light- like Wilson lines U(x±) are formally Mobius invariant and this 
is why the leading-order BK equation is conformal. However, because our cutoff of the rapidity divergence is not 
invariant, the NLO evolution kernel (|44p has the non-invariant double-log term. To illustrate the non-invariance of 
the dipole with the cutoff (jHJ) let us consider the high-energy operator expansion Q of two BPS-protected currents 

O = Yf* Tl '{ z2 } (Z = -7ni<t>i + ifa)) in color dipoles. The Regge limit of the amplitude 

A(x,y;x',y') = (x - y)\x' - y'f{0 (x)0(y)0(x')0(y')) (80) 
can be achieved by the rescaling Q 

x -> \x*p! + -x,p 2 +x±, y -> Ay* + jV*P2 + V± (81) 

with A — > oo. In this regime, the T-product of the currents 0(x) and 0{y) can be expanded in color dipoles as follows: 

T{6{x)6{y)} = Jd 2 z 1 d 2 z 2 I h °(z 1 ,z 2 )Tr{UlU^} 

+ J d 2 Zl d 2 z 2 d 2 z 3 I Nh0 (z 1: z 2 , z s )[Tr{T n UlU%T n UlU%} - N c Tr{UlU%}] (82) 

(structure of the NLO contribution is clear from the topology of diagrams in the shock-wave background, see Fig. [7| 
below) . 

Let us calculate the impact factor taking x, = y. = for simplicity. The leading-order impact factor is proportional 
to the product of two propagators (|55|) 
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where Zi 



and the color trace is taken in the adjoint representation. This expressions coincides 



with the result of Ref. [22j- It is easy to see that under the inversion x± — > x_l/xj_,x* — > x*jx\ 



x± -> xi_/x\,x* -> x*/x 2 ± , yj_ -> y±/x 2 ± ,y* -> y*/yj 



(84) 



^ is transformed as -E^ — > 2^ 4 iJi so the leading-order impact factor is invariant. 

The NLO impact factor for two Z 2 currents is given by the two diagrams shown in Fig. [Jj To calculate them we 





(a) (b) 

FIG. 7: Diagrams for the NLO impact factor, 
will use the following representation of scalar and gluon propagators in the shock-wave background (at > > ) 
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where 4- can be either » \ . or „ . which leads to the same result after subtraction of the leading-order contribution 
(see below). 

To calculate the next-to-leading impact factor we need the three-point scalar-scalar-gluon vertex Green function 
(vertex with tails): 
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Since zi* — z» = one can easily check that #(z») in the l.h.s. of this integral can be erased since the contribution 
from z* < vanishes. After that, one easily gets Eq. (|87| from the conformal integral 
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which is easily calculated by inversion z^ —> z^j z 2 
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Now we are in position to calculate the contribution of the diagram in Fig. [7Ji: 
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Using the Eq. (|87j) one can reduce this equation to 
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Performing integration with respect to z.j we get 



2n 2 [(x-z 3 ) 2 -ie]d 3 , 27r 2 [(y-z 3 ) 2 



(T{6(x)6(y)}f^ 



2a s f d 2 z\d 2 z 2 

(N 2 -l)^xly 2 J Z 2 Z 2 



d 2 z 3 



(Z13,Z 23 )^ 



■Tr{T n U Zl UlT n U Z3 U}j 

Jo 



" da 



(recall that Z\ 



(T{6(x)6(y)}f2 M = 



^13 ^23 

Similarly, one can demonstrate that the diagram shown in Fig. [7t> yields 
2a, 



d 2 z\d 2 z 2 



1^-2 J z 13 JO a 



(N 2 - l)^x 2 y 2 

Using the z\ «-> z 2 symmetry in this equation one can write down the sum of Eq. (|91[) and Eq. (|92[) as 



f^3 



(T{6(x)d( y )}f* m 



d 2 z\d 2 z 2 



(iV 2 - l)7r% 2 y 2 
It is convenient to rewrite this in the following form 



-72-72 



d 2 z 3 



z 12 
Z 13 Z 23 



Tr{T n U Zl UlT n U Z3 Ul 2 } / 
Jo 



' da 



giaf Z 3 



(90) 



(91) 



(92) 



(93) 



(T{G(x)0(y)}) A 
- N c Tv{U Zl Ul} 



FigE] = 

da 



ia$Z 3 



1)~ x *y*J Z/ l- & 2 J Z 13 Z 23 

^1^2 J z 13 z 23 JO 



2a s N c 



(N 2 - l)n*xly 



'da 



ia^Z 3 



(94) 



Let us discuss now the contribution of Fig. [5] diagrams. Since this contribution is proportional to Tr{U Zl Ul } it 





(a) 



(b) 



FIG. 8: Diagrams for the NLO impact factor without gluon-shockwave intersection. 

can be restored from the comparison of Eq. (|94|) with the pure perturbative series for the correlator (T{0(x)0(y)}). 
Indeed, if we switch off the shock wave the contribution of the Fig. [7] diagrams is given by the second term in Eq. 

(with U,U^ repaced by 1). On the other hand, perturbative series for the correlator (T{0(x)0(y)}) vanishes 
|23j | and therefore the contribution of the Fig. [5] diagrams should be equal to the second term in the r.h.s. Eq. ([M]) 
with opposite sign. Thus, the first term in the r.h.s. Eq. (|94|) gives the total contribution to the impact factor: 

(T{6(x)6(y)}) F A ism ^ Fis ® 

The integral over a in the r.h.s. of Eq. (|95p diverges. This divergence reflects the fact that the r.h.s. of Eq. (j9"o]) 
is not exactly the NLO impact factor since we must subtract the matrix element of the leading-order contribution. 
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Indeed, the NLO impact factor is a coefficient function defined according to Eq. ([82]) . To find the NLO impact factor, 
we consider the operator equation (|82|) in the shock- wave background (in the leading order (U Z3 )a = U Z3 ): 



(T{0(x)0(y)}) A - \d 2 z x d 2 z 2 I LU (x,y;z 1 ,z 2 )(Tr{UlU^}) 



= J d 2 z 1 d 2 z 2 d 2 z 3 I NLO (x, y; zi, z 2 ,z 3 ; r,)[Tv{T n U Zl UlT n U Z3 Ut 2 } - N c Tr{U Zl Ul 2 }] 
The NLO matrix element (T{d(x)G(y)}) A is given by Eq. $Q while 
'd 2 Zl d 2 z 2 I LO (x, y; z u z 2 )(Tr{UlU%}) A 



(96) 



7T 2 (iV 2 



_ I d 2 Zld 2 Z2 f^Lf d 2 Z3 J^ [TT{T "U Zl UiT n U Z3 U! 2 } - N c Tr{U Zl Ul}} 

L)J X, 1 X, 2 TT J a J Z \ 3 Z 23 



(97) 



as follows from Eq. (116p . The a integration is cut above by a — e 71 in accordance with the definition of operators U v 
©. Subtracting flST]) from Eq. ([55]) we get 



rNLO 



(x,y;z 1 ,Z2,z 3 ;T]) 



a s (x*y*Y 



- 13 



-13 



^(N 2 -l)z 2 2 z 2 3 Z 2 Z 2 



, as ITT ■ 

In — Z 3 YC 

4 3 2 . 



! da 



a 



da 
a 



(98) 



Let us rewrite the operator expansion (|82|) in the explicit form 

(x-y) 4 / j2 _ , 2/v {x*y*)~ 2 



{x-yfT{0{x)0{y)} 



tt 2 (N 2 1) 



d z 1± d z 2 ±- 



Z 2 Z 2 



(99) 



a s (x - y) 4 
n\N 2 - 1) 



d zi±d z 2 ±d z 3 



Z] 2 (•£*£/* ) 
z 13 z 23^1^2 



, as ITT ■ 

Xn-Z 3 -- + C_ 



It is easy to see now that under the inversion |84J the leading-order impact factor is invariant while the NLO impact 
factor is not because of the non- invariant logarithmic term ln^iS^. Since the original T-product of the currents in 
the l.h.s. of Eq. (|9"9")l is conformal, it indicates that our operators U v with the "rigid cutoff" © are not Mobius 
invariant. However, if we expand the original T-product in composite conformal operators ([7]) instead, the resulting 
impact factor is conformally invariant: 

( X ~ V) i .,2 . ,2. •":>• r r ( , .„ , , |,.,.„i ( | 00 ) 



(x - y) 4 T{0(x)0(y)} = 



^{Nl - 1) 



d 2 z 1± d 2 z 2± { -^lTl [rrifyug}]"** 
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^\'°2 



a s {x-y) 4 
2^{N 2 _ r 

2^2 J2 



d z 1± d z 2 ±d z 3 



z i 2 { x *y*) 

r 2 ..2 7-2 72 



■'13' i 23^1 ^2 



x In 



a<7 s z± 2 



16z 2 3 z 2 3 



(a - z 3 f (y - z 3 f 



ztt + 2C) [Tr{T n UlU%T n UlU%] - N c Tr{UZU%}] 



The arbitrary dimensional constant a should be chosen in such a way that the impact factor in the r.h.s of Eq. (1101|) 
does not change under the rescaling (f8"Tj) . The proper choice for our T-product is a = s2 ^^ so our final operator 



expansion takes the form 

(x - y) 4 T{6(x)6(y)} 
Oi s {x - y) 4 



{x - y) 4 
* 2 {N 2 c -l) 



d Z\d z 2 



jx*y*) 2 

7 2 7 2 



conf 



(101) 



^ 4 (N 2 - 1) 



d z\d z 2 d z 3 



z i 2 (x*y*) 



hi 



2 p 2 V 



I6{x - y) 2 j_z 2 3 z 2 3 



(x - z 3 ) 2 (y - z 3 ) 2 



x [Tr{T n UlU%T n UlU%} - N c Tr{UlU%}] 



where the conformal dipole [T^{U Zl Ut 2 }] coni is given by Eq. ([7]) with a — . Now it is evident that the impact 

factor in the r.h.s. of this equation is Mobius invariant and does not scale with A so Eq. (|4"5|) gives conformally 
invariant operator up to a 2 order. In higher orders, one should expect the correction terms with more Wilson lines. 
This procedure of finding the dipole with conformally regularized rapidity divergence is analogous to the construction 
of the composite renormalized local operator by adding the appropriate counterterms order by order in perturbation 
theory. 
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C. Leading-order evolution of the four- Wilson-line operator 

In this Appendix we derive the evolution equation for the four- Wilson-line operator r Tr{T a U^U^T a U^ 3 U^} in the 
leading order in perturbation theory. As a first step, we rewrite the hierarchy equations of Ref. [5j in the adjoint 
representation. It takes the form 

i-0l)iAUl)u = S / rf2 *3 %^[(T ^J«(^)« + (^T 6 )«(r ^)K](2^-^-^)- 6 ) (102) 

' J 13 23 

for the pair-wise interaction and 



, <Kh = "i / —(T a UlT%(Ul - Ul) ab (103) 



for the self-interaction. Using the above equations we get 



^ [ d 2 z 4 (UD aa ' \ ^ Z ^ Tr{T a T b Ul [T a ' , T b '}U^ + [T b , T a \U\T b ' T a ' U^iW^ - XJ\ - U2J bb ' 
2ir J L z 14 z 34 



7 2 2 
z 24 z 34 

(Zl4, z 2i) 



~r~^—Tr{[T a ,T b }UlT a 'T b 'u^ + T b T a Ul[T b ' ,T a ']U^}(2Ul - IJ\ - Ul d ) bb ' 



-14^24 



Tr{T a T b U!>T a 'T b 'u^ + T b T a U^T b ' T a 'u^}(2U^ - fj\ - U^) bb ' 



- ^Tr{{T a ,T b ]Ul[T a \T b '}U^}(2UI> 4 - 2U2J bb ' 

Z 34 

+ -^TT{T a T b U r JT b 'T a 'u^}(2Ul - 2UI> 1 ) bb ' + l-TY{T b T a U^T a ' ' T v 'U^}(2U^ - 2U2 2 ) bb ' 

Z 14 ' ' Z 24 

m AT f 7 2 

U Zi ^{TT{T«U Zl UlT«U z JJl} - N c Tr{U Zl Ul 2 }) (104) 

" J z 14 z 24 

After some algebra, the r.h.s. of Eq. (] 1 04[) reduces to 



Ota 
4^ 

z 13 
7 2 2 
z 14 z 34 

z 23 
7 2 2 
z 24 z 34 

a s 7V c 



d 2 Zi (U2 3 ) aa ' -^ r Tr{T a T b U!! 1 T a 'T b 'u}jJ + T b T a U2 1 T b 'T a 'u^}(2Ul ~ - ^ 2 ) b6 ' 

Lz 14 Z 24 

Tr{T Q T b L> z r ' i [T a ',T b ']L> z t 2 " + [T b ,T a }UlT b ' T a ' U^}(2U r J 4 - U\ - Ul{) bb ' 
Tr{[T a , T b \U\T a ' T b ' U\2 + T b T a U\ [T b ' ,T a ']U^}(2Ul - &1 ~ Ul) bb ' 
d 2 z,^^{T?{T a UlUt2T a UlUi?} - N c Tr{U Zl Ul 2 }) (105) 

z 14 z 24 
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D. Color traces 

In this section we rewrite the adjoint traces in our evolution equation ([6]) in terms of the fundamental traces. The 
master formula for traces has the form (in this section we use the space-saving notations Ui = U Zi ) 

Tr{[T a , T b ]U?T a V ! Xj\ + T b T a U\ [T b> ,T a ']UljUg a 'uf (106) 

= \ [ - tTiUiU^tT^UljtiiUiUljtriUsUl} - ^{UlU^tviUlU^tviulu^tviUtUi} 

+ tr{UiUl}tT{U 2 UlU 3 ulUiUl} + tv{UxUl}tx{U 3 ulu 4 UlU 2 ul} + tx{U 3 ul}tx{U 4 U 2 j UiU%U 2 Ul} 
+ ix{U2Ul}tr{U4.U}U 3 UlUiU£} - tr{UiUl}ix{U4.U^U 3 UlU 2 Ul} - tr { C/ 2 ^ } tr { t/ 3 ^ C/ 4 t/ 2 f C^i C/l } 

- triU^UiUljtriul^UlUs} + tx{UlU 2 }tx{UlU A UlUxUlU 3 ) + tr{ C/^ C/ 3 }tr{ £/ 2 C/2 C/i C/ 4 } 
+ tr^l^M^U^CWl^i} + ixiUlUtytiiUlUtUlfyUtUa} - tr{{/ 1 t {/4}tr{[/ 4 t [/ 2 {/ 3 t {/iC/ 2 t [/3} 

- tr{U^U 3 }tr{U^UiUlU 2 UlU 4 } - tr { C// 1/ 2 1/] C/ 3 } tr{ i7i f7+ i7 4 C/J } 

where Tr and tr stand for the trace in the adjoint and the fundamental representation, respectively. We will also need 
traces 

Tx{T a UiT b Ul}Uf = Tr{T a J7i[/ 3 t r b C/ 3 C/ 2 t } 

= \te[{UiUl}tT{U 2 Ul}tv{U 3 Ul} + ix{UiU£}tt{U a Ul}ix{U2Ul} - ix{UiU%UaUlU 2 u£} - triUiU&UlUsU}}], 

Tr{[T Q ,T b }U 1 T a ' ' T v ' u\ + T b T a Ui[T b \T a ']ul}U% a> 'uf 
= Tr{[T a , T b \UlT a ' T b ' U\ + T b T a U^[T b ' ,T a ']ul}Uf 'uf 

= Tr{[T a ,T b ]U , 1 l T a 'T b 'ul +T b T a U^[T b ' ,T a ']U%}Ug a ' = - N c Tr{T a U 1 UlT a U 3 Ul}, 
Tr{T a T b UiT a ' ' T b ' ' Xj\ + T b T a UiT b ' T a ' ut,}U? a ' l/f 

= Tr{T a T b U 1 T a 'T b 'u^ +T b T a U 1 T b 'T a 'u^}U^ a 'ul b ' = N c Tx{T a UiUlT a U 3 Ul} (107) 

which easily follow from Eq. (|106[) . 

Let us now find the master trace (|106[) in the two-gluon (BFKL) approximation. First, note that in this approxi- 
mation the trace of six Wilson lines reduces to one dipolc. Indeed, 

tr{UiUlu 4 U$U 3 ul} = tr{(f7iZ7 3 f - 1)U 4 (U%U 3 - l)u\} + tr{(lW 3 f - 1} + tx{(U$U 3 - 1} + 

= tr{(^ 1 ^ 3 t -l)^ 3 (C/ 2 t ^ 3 -l)^ 3 t } + tr{(C/iC/ 3 t -l} + tr{(J7 2 t ?7 3 -l} + 7V c 2 = tx{UiU%} (108) 

where we replaced by t/4 by U 3 in the first term in the second line since it does not matter in the two-gluon approx- 
imation (both of them should be replaced by 1). All other six- Wilson-line terms in the r.h.s. of Eq. (| 1 06[) can be 
similarly reduced to single-dipole terms. 

Using the same trick we can reduce the trace of four Wilson lines to the sum of dipolcs: 

tv{UxUlUiUl} = tr{(Lr 1 U 2 t -l)(^ 4 ?73 t -l)} + tr{?7iU 2 t -l)} + tr{U4C4-l)} + ^c ( 109 ) 
= txiiUxUl - l)U 2 U\{UiUl - 1)} + tv{U x Ul - 1)} + tr{U 4 ul - 1)} + iV 2 

= tx{{U x Ul - 1)} + tv{U 2 Ul ~ 1} + ti-{U A ul - 1)} + trjLWj - 1)} - tx{U x Ul - 1)} - tx{U 2 U% - 1)} + iV 2 

where in the first term in the second line we have inserted U 2 U\ since it does not matter in the two-gluon approxi- 
mation. 

Using formulas (|108[) and (|109p it is easy to demonstrate that all the terms containing traces of four and six Wilson 
lines in the r.h.s. of the Eq. (|106[) sum to 2iV 2 in the two-gluon approximation so one gets 

Tr{[T a , T b ]U^T a ' T b ' Xj\ + T b T a U^ [T b \T a ' pljuf uf 2 ~ 9 = °" - A/ 2 (iV 2 - 1) - ^iV c 3 [tr^Kf - 1} 
+ tx{u\u 2 - 1} + tx{U 2 u\ - 1} + tx{UlUi - 1} + tx{U A Ul - 1} + tx{U\U 3 - 1} + tx{U 3 Ul - 1} + tx{U%Ui}] 
= - iV 2 (iV 2 - 1) + \n 2 c {N 2 c - 1) [V(z u z 2 ) + V(z u z 3 ) + V(z 2 , Zi ) + V(z 3 , z 4 )] (110) 
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where we use the definition ((54)) . For completeness, let us present also Tx{T a UiUlT h U 3 Ul} in the two-gluon approx- 
imation. From the first line in Eq. (|10T|) one easily obtains 

^{T a U x UtT b U 3 Ut} 2 ~ 9 =° n N C (N 2 - 1) - \n c (N 2 c - 1) [V( Zl , z 2 ) + V(z u z 3 ) + V(z 2 , z 3 )] (111) 
To describe the conformal operator (|48p we need one more trace (which eventually drops out after the integration of 

2 

the operator (j48|) over z 3 with the weight 2 12 a ) 

Z 13 Z 23 

Tr{T a T b U 1 T a 'T b 'u^}U^ a 'uf = - tx{U x Ul}tx{U 3 U ] 2 }tx{U 2 Ul}tx{UiUl} (112) 

+ ixiUsUhtiiUiU&triUiU&tiiUlU!} + ^{UzUlUiU^triUtUlUzUl} + tr{ (7 2 ^ i7i t/] } tr { C// C/ 3 C/l ^4 } 
+ triUiU^triUsUl^ulUiU}} + ix{UiU^}ix{UlUaUtUil4u 2 } + tr^^M^I ^4^1*4} 
+ tr^^M^l^i^^} - tr{^ 3 V{ W^^t^l} - tr^^M^^^M^} 

- tr{l/3l^}tr{I7^2l/]^it^C^} - tr^Ul^riUxUlU^UlUiUl) - tr^^ML^^M^} 

- tx{U A Ul}tx{UlU 3 ulu x ulU2} - tx^U^txiUxUlUsUlUiUl} - te{L^L7i}tr{U^L^Z7 3 L/J} 

In the two-gluon approximation this yields 

= \n 2 c (N 2 c - 1) - \n 2 (N 2 - 1) [V(z u z 3 ) + V(z 2 , z 3 ) + V{z uZi ) + V(z 2 , z 4 )] (113) 



E. Integrals 



In this section we calculate two basic integrals which we use in this paper. The first one is 



,2 r d 2 Z4 r ,2 .2 ,2 , ,2 .2 



-12 



13^23 



'13 



-23 



-12 



'14*34 



'24^34 



'14*24 



ln^ft = 4C(3)[5(«a3) - S(z 13 )} 



'14*23 



The easiest way to prove this at z 3 ^ Z \ , z 2 is to set z 2 — and make an inversion x 
reduces to 



(114) 

1/x so the integral (|114p 



2(£i-z 2 r d 2 hji 



(h - Z 3 , Zl 



*) IngL^l! = o 



(Si - z A ) 2 (z 3 - z 4 ) 2 (zi - z 4 ) 2 
The 5-function terms in the r.h.s. of Eq. (|114p can be restored from the formula 



d 2 z 3 d 2 Zi 



(Z12,Z23) ( 1 



-31 



y z 13 z 24 



"23*14 



^12 (Z12,Z23) 



-13*24^14*23 



hi 



-13*24 



'23*14 



4C(3) 



which follows from the integral 



,d d 



(x,z) 



(x,z) 



x 2 (x, z) 



(x- zf(z- z') 2 z' 2 (x- z') 2 (z- z'Yz 2 {x- z) 2 z 2 (x~ z') 2 z l2 \ {x~z') 2 z 



In 



(x — z) z 



2^/2 



v 2' 2 ' v '2 ' (x 2 ) 2 - d X ^ v 2 



- 2<ip(d - 2) + 2^(1) - 2^(2 - -) 



*<|) 



r(if)r(3^-2)r 2 (2-C 



d^2 



d-2 r 2 (d- i)r(3 - rf) J 

The second integral is somewhat more tedious 



4C(3) 



-12 



d 2 Z4 



Z l 3 Z 2 3 J " 



-23 



In 



222 

z 12*13*34 



6 13 



In- 



'34*24 



-14*24 



-14*34 



-13*23 



'12*14*34 



-12 



'14*24 



In 



-14*24 



'13*23 



= 2 



_fl2_ l n fl 2 . ln fl2 

r 2 r 2 111 -2 111 r 2 

13 23 •-'i^ **o 



8C(3)<5(z 13 ) 



6 13 



-23 



(115) 



(116) 
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First we prove this equation at z% ^ Z\,z<i. To simplify the notations, we take zi = an denote z\ = x, Z3 = z and 
z 4 = z'. Let us start from the first term in the square brackets in the l.h.s. of Eq. Q116P 



d 2 z' 



■In 



x 2 (x — z) 2 (z — z') 2 



{z-z'fz' 2 ix-z'Yz' 2 



(117) 



2 / dV 



(z - z') 2 z' 2 



111 



(x — z') 2 , .<- 



In- 



2 . / 



cTz 



{x-z) 2 J ' (z-z') 2 z' 2 J " (z-z') 2 z' 2 X 



1 (z,z-z') , (x-z') 



4 / d 2 z' H^ln- 



A2 



We get 



-2 / dV 



(z - z>) 



= lim I - 2 / d d z' 



M2 „/2 



hi 



- z' f 



-In- 



d 2 ~' 



(x - z) 2 J (z — z') 2 z 



2 r /2 



(118) 



1 



1 



- *'12 ,/2 



(z - z') 



In 



(a: - z') 



A2 



-In- 



(x~z) 2 J (z-z') 2 z' 2 



lim 27r 

d^2 



r(i-|) r(i-4) r(2-|) 



In 



and 



((x-z) 2 ) 1 -! (x 2 ) 1 ^ (z 2 ) 1 -^ (x-z) 2 
/ {z,z-z') , (x-z 1 ) 2 



B 



d d 
2' 2 



-4 / d d z' In 

(z- z') 2 z' 2 

- - I j In :- +-(- - I j ( In 



TT d ' 2 T(d/2) ,d , , d x 

lim 4 / 7 - 1 r 2 - - 

m d-2rd-l l 2 ' y 2' 



, (a; — z) 2 x 2 , (x - z) 2 
7r In In - 



1 lnz 2 



2 V2 

2 7 rlniln ( ^ z)2 



-In 2 



ln 2 x 2 ) -1- ( - - 1 ) ln^ - - ( - - 1 1 In 



(x — z) 2 (x — z) 



1 id 



2\2 



2 2 



so 



d 2 z' 



z 2 x 2 (x — z) 5 



hi 



(z-z') 2 z' 2 (x-z') 



7rln" 



(x - z) 2 



(119) 



(120) 



Similarly 

d 2 z' 



(x — z) 2 x 2 (x — z') 2 (z — z 1 ) 2 



In 



(x — z') 2 (z — z') 2 [x — z)^z 



■ hi 



(x — z') z 



2„/ 2 



lim j?r r(2 i)r 2 ( d 2 ) c _ d _ 2 in (x zfz 2 _ _ d _ 2 

1 L 



(x-z') 2 z' 2 (x-z) 2 z 2 

1 



, . ... , _,. + ul (x-z) 2 -0(2--)+0(-) 

d^2 (| - l)r(d- 1) V a; 2 L d-2 V ' n 2> VK 2> 

+ 0(1) - i>{d - 1)] + 2|x| d - 2 [ - + lnx 2 - 0(2 - ~) + ^(|) + 0(1) - 0(d - 1)] - |x| d - 2 ln(x - z) 2 z 2 } 



2 2 2 

i2 rt l l *^ 

7T In — + 2n In — In — 

(x — z) z z z [x — zy 



(121) 



and therefore 



2J 



<Tz 



hi 



(x - z') 



A4 



(x-z) 2 



In- 



(x — z) 4 z 2 



, (x — z') z 



2,/ 2 



In 



(z-z') 2 z' 2 x 2 (x-z) 2 (x-z') 2 (z-z') 2 x 2 (x-z') 2 (z-z') 2 (x-z') 2 z' 2 (x-z) 2 



= 27rln^5-ln- 

z z (x — z) z 



(122) 



However, it is easy to see that some sort of 5-function contribution to the r.h.s. is necessary. If we integrate the l.h.s. 
over z with the weight ^ x ^ z yi z i we get zero because of the antisymmetry of the integrand with respect to z z'. On 
the other hand, the integral of the r.h.s. does not vanish because 



d 2 z 



In — In ■ 



(x — z) 2 z 2 z 2 (x — z) 



4ttC(3) 



(123) 
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To fix the coefficients in front of possible 5-function contributions ~ S(z) and/or S(x — z) we calculate the integral of 
the l.h.s. of Eq. (fT22l with the trial function ■ We get 



— ^ / d 2 zd 2 z 



In 



z 2 (x — z) 2 



■In- 



z A (x — zY 



z 2 (x~ z) 2 l{x~ z') 2 (z~ z') 2 [x~z') 2 z' 2 z' 2 {x-z') 2 {z-z') 2 z' 2 x 2 z' 2 {z-z') 2 



lim B(i, t - l)B(d -1,1- + H~ 2 ) + - d) + ^1 - 2) - 2*(d- 1) - 2^(2 - |)} 



r(3 - d) 



d^2 x 2 2 
4C(3) 



(124) 



where B(a,b) — T(a)T(b)/T(a + b). It is clear now that the result for the integral in the l.h.s. of the formula (|116|) 
should be as cited in the r.h.s. of Eq. (|1 16[) - it satisfies both the Eq. (|124[) and the requirement that the integral of 



the l.h.s. with the trial function 



vanishes. 



z 1 {x — z)' 1 

We will need one more integral which is obtained by antisymmetrization of Eq. (I116P with respect to Z\ <-> z^ 



'12 



'13^23 



d 2 Z 4 



'13 



In 



-12 ^34 



-13^-24 



2 2 2 

Z 23 i Z 12 2 34 
> n m o n 



'14^-23 



= 4C(3)[<S(Z13)-*(Z23)] 



(125) 
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